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Abstract 

In this paper we propose a simple and efficient method to compute the ordered 
default time distributions in both the homogeneous case and the two-group het- 
erogeneous case under the interacting intensity default contagion model. We give 
the analytical expressions for the ordered default time distributions with recursive 
formulas for the coefficients, which makes the calculation fast and efficient in finding 
rates of basket CDSs. In the homogeneous case, we explore the ordered default 
time in limiting case and further include the exponential decay and the multistate 
stochastic intensity process. The numerical study indicates that, in the valuation 
of the swap rates and their sensitivities with respect to underlying parameters, our 
proposed model outperforms the Monte Carlo method. 
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1 Introduction 



Modeling portfolio default risk is a key topic in credit risk management. It has important 
applications and implications in pricing and hedging credit derivatives as well as risk 
measurement and management of credit portfolios. There are two strands of literature 
on credit risk analysis, namely, the structural firm value approach pioneered by Black 
and Scholes (1973) and Merton (1974), and the reduced-form intensity-based approach 
introduced by Jarrow and Turnbull (1995) and Madan and Unal (1998). In the classical 
firm value approach the asset value of a firm is described by a geometric Brownian motion 
and the default is triggered when the asset value falls below a given default barrier level. 
In the reduced-form intensity-based approach, defaults are modeled as exogenous events 
and their arrivals are described by using random point processes. 

The reduced-form intensity-based approach has been widely adopted for modeling 
portfolio default risk. Two major types of reduced-form intensity-based models for de- 
scribing dependent defaults are bottom-up models and top-down models. Bottom-up 
models focus on modeling default intensities of individual reference entities and their 
aggregation to form a portfolio default intensity. Some works on the bottom-up ap- 
proach for portfolio credit risk include Duffie and Garleanu (2001), Jarrow and Yu (2001), 
Schonbucher and Schubert (2001), Giesecke and Goldberg (2004), Duffie, Saita and Wang 
(2006) and Yu (2007) etc. These works differ mainly in specifying default intensities of 
individual entities and their portfolio aggregation. Top-down models concern modeling 
default at portfolio level. A default intensity for the whole portfolio is modeled without 
reference to the identities of the individual entities. Some procedures such as random 
thinning can be used to recover the default intensities of the individual entities. Some 
works on top-down models include Davis and Lo (2001), Giesecke and Goldberg (2005), 
Brigo, Pallavicini and Torresetti (2006), Longstaff and Rajan (2007) and Cont and Minca 
(2008). 

One of the major applications of portfolio default risk models is the valuation of credit 
derivatives written on portfolios of reference entities. Typical examples are collateralized 
debt obligations (CDOs) and basket credit default swaps (CDSs). The key to valuing 
these derivatives is to know the portfolio loss distribution function. The kth to default 
basket Credit Default Swap (CDS) is a popular type of multi-name credit derivatives. 
The protection buyer of a kth to default basket CDS contract pays periodic premiums to 
the protection seller of the contract according to some pre-determined swap rates until the 
occurrence of the kth default in a reference pool. Whereas, the protection seller of the kth 
to default basket CDS pays to the protection buyer of the contract the amount of loss due 
to the kth default in the pool when it occurs. Different approaches have been proposed in 
the literature for pricing the kth to default basket CDS under the intensity-based default 
contagion model. Herbertsson & Rootze (2006) introduce a matrix-analytic approach 
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to value the kth to default basket CDS. They transform the interacting intensity default 
process to a Markov jump process which represents the default status in the portfolio. This 
makes it possible to apply the matrix-analytic method to derive a closed-form expression 
for the kth to default CDS. Yu (2007) adopts the total hazard construction method by 
Norros (1987) and Shaked & Shanthikumar (1987) to generate default times with a broad 
class of correlation structure. He also compares this approach with the standard reduced- 
form models and alternative approaches involving copulas. Zheng & Jiang (2009) use the 
total hazard construction to derive the joint distribution of default times. They give a 
closed-form expression for the joint distribution of the general interacting intensity default 
process and an analytical formula for valuing a basket CDS in a homogeneous case. 

Here we propose a simple and efficient method to derive the kth default time distri- 
bution under the interacting intensity default contagion model. We give the recursive 
formulas for the ordered default time distributions, and further derive the analytic solu- 
tions in a group of homogeneous entities and in two groups of heterogeneous entities. In 
the homogeneous case, we discuss the ordered default time in limiting case and further 
include the exponential decay and the stochastic intensity process. We derive the pricing 
formula under a two-state, Markovian regime-switching stochastic intensity model. In 
addition, we show that our proposed method is superior to the simulation method in 
studying the sensitivities of the swap rates to changes of underlying parameters. 

The rest of the paper is organized as follows. Section 2 gives a snapshot for the in- 
teracting intensity-based default model. Section 3 discusses the homogeneous case and 
applies the recursive method to characterize the ordered default time distributions de- 
rived in Zheng & Jiang (2009). Section 4 extends the method to study the multi-state 
stochastic intensity process. Section 5 addresses the two-group heterogeneous case. Sec- 
tion 6 presents numerical experiments for the evaluation of the basket CDS under various 
scenarios and the sensitivity analysis. Section 7 concludes the paper. 

2 A Snapshot for Interacting Intensity-Based Default 
Model 

In this section, we give some preliminaries for the paper to facilitate our discussion. Let 
(Q, J 7 , {J-i}f>o,-P) be a complete filtered probability space, where we assume P is a risk- 
neutral martingale measure (such a P exists if we preclude the arbitrage opportunities), 
and {J-t}t>o is a filtration satisfying the usual conditions, (i.e., the right-continuity and 
P-completeness). We consider a portfolio with n credit entities. For each % — 1, 2, • • • ,n, 
let n be the default time of name i. Write Ni(t) = l{ Ti <t} for a single jump process 
associated with the default time n, and {J~t~}t>o is the right-continuous, P-completed, 
natural filtration generated by iVj. Suppose {X t } t > is the state process, which represents 
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the common factor process for joint defaults. Write {J-*}^ f° r the right-continuous, 
P-completed, natural filtration generated by the process {X t } t > . For each t > 0, write 

jr t = jf v 7} V . . . V . 

Here Tt represents the minimal a-algebra containing information about the processes X 
and {iVj}™ =1 up to and including time t. 

We assume that for each i — 1, 2, • • • , n, JVj possesses a nonnegative, {J-"t}t> -predictable, 
intensity process Aj satisfying 

E ^ Xi(s)ds) < oo, t > 0, 
such that the compensated process: 

Miit) := Ni(t) - / X^ds , t > , 
J o 

is an ({J-i}t>o, P)-martingale. We further assume that the stochastic process X is "ex- 
ogenous" in the sense that conditional on the the whole path of X, (i.e., J 7 ^), Aj are 
{Vi J~\ }t>o-predictable. 

To model the interacting intensity default process, we consider the following form: 

Xi(t) = Oi{t) +Y,bij(t)e~ dii{t ~ Tj) hT j <t}, i = 1,2, . . .,n, (1) 

where and %(£) are J-" x -adapted processes, and dij are positive constants represent- 
ing the rates of decay. The introduction of exponential decay into the intensity-based 
default process is of practical significance, which indicates once a default occur in the 
portfolio, its effect on the other surviving entities will decrease at a rate proportional to 
its present impact. Intensity rate processes in (1) determine the probability laws of the 
default times. To price the kth to default basket CDS, the distribution of the kth default 
time has to be known. 



3 Homogeneous Case 

In this section, we present a simple method to derive the distribution of the kth default 
time in a group of homogeneous entities under the interacting intensity default model. 
Our method is based on the kth default rate and the distribution of the random duration 
between two defaults, where the contagion intensity process follows: 

Hi) = a^l + £ ce- d ^h {Tj < t} j , i = 1, 2, . . . , n , (2) 

where a is a positive constant and c and d are nonnegative constants. 
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We note that the Markov chain approach (Herbertsson & Rootze(2006)) cannot solve 
this kind of processes with exponential decay (d > 0). Zheng & Jiang(2009) adopt the 
total hazard construction method to give the joint distribution of default time Tk, k = 
1,2, ... ,n, while finding the ordered default time distributions remains to be intractable. 
Here we give the recursive formula of the joint distribution of the ordered default times 
by our proposed method. Let X h+1 (t) be the (k + l)th default rate at time t that triggers 
T h+1 . Then X k+1 (t) will be the sum of the default rates of the surviving entities after r k . 
Under the homogeneous situation, given the realization of r l , i — 1, . . . , k, the {k + l)th 
default rate is: 

k 

\ k+l (t) = a(l + Y,ce- d(t - Tl) )(n - k), 
i=i 

where r k <t < r k+1 . Then we have 

P(r k+1 - r k > t | r\ i = 1, . . . , k) = exp{- / X k+1 (r k + s)ds} 

Jo 

Thus, 

P(r k+1 > 1 1 t\ i = 1, . . . , k) = exp{- C X k+1 (s)ds}, 

Jr k 

where t > r k , which implies 

U^,i=i,...,k{t) = a{l+j^ce- d ^){n-k)e W {- f a{l+j^ce- d ^){n-k)ds], (3) 

i=i Jrk i=i 

where f T i (t) = nae~ nat . One can apply (3) to derive the joint density function of r 1 , r 2 , . . . , r k+1 
with the following recursion: 

•/r 1 ,T 2 ,...,T fc + 1 (^l) h, ■ ■ ■ , tk+l) = /r fe + 1 |r i =t i ,i=l,...,fe(^fc+l)/r 1 ,r 2 ,...,T fe (^l; h, ■ ■ ■ , tk) , 

where ti < t 2 < ■ ■ ■ < tk+i- The unconditional density function of r fc is given by the 
integral: 

ft rtk-i rt-2 

f T k(t) = / / •••/ f T i r 2 ,..., r fc-i ,r k (ti, . . . , tft-i, t)dt\ ■ ■ ■ dtk-2dtk-i (4) 
Jo Jo Jo 

Example 1 If n = 2, then the joint density function of r 1 and r 2 is: 



f T l T 2{tl,t 



2) 



2a 2 (l + ce-^ 2 -^) exp{-a(t! + t 2 ) + ^(e"^-* 1 ) - 1)}, t x < t 2 
0, t x > t 2 . 



The unconditional density function of r 1 is given by 

f Tl (t) = 2ae- 2at , 

and that of r 2 is given by 

fAt) = 2a 2 f\l + ce^*-* 1 )) exp{-a(t + t t ) + ^(e"^-^ - \)}dt x . 
Jo d 
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If we simplify our model by assuming that d — 0, we have the following proposition. 

Proposition 1 Suppose there are n entities in our portfolio, where the contagion intensity 
process follows 

\i{t) = a ll + J2 cl {Tj<t}\ , z = l,2,...,n. (5) 

V w J 

The kth default time r k is the sum of k independent exponential random variables, i.e., 

k-l 
i=0 

where Xi are independent exponential random variables with rates a(l + ic)(n — i) respec- 
tively. The the unconditional density function of r k , k=l,2,. . . ,n, satisfy the following 
recursive formula: 

f T k+i(t) = ail + kc)(n - k) f f Tk {u)e- a(1+kc){n ~ k){t ~ u) du, (6) 

Jo 

where the initial condition is: 

f T i{t) = nae~ nat . 

Proof: We note that r k+1 — r k is independent of t 1 for % = 0, 1, . . . , k, where r° is 
assigned to be 0. Due to the homogeneous and symmetric properties of the entities in our 
portfolio, we have 

\ k +\t) = a{l + kc)(n - k) for r k < t < r k+1 . 

Therefore 

p( T k+1 — r k > t) = e - a ( 1+fcc )(«- fc )* 

and this implies 

f T k+i_ Tk (t) = a(l + kc)(n - k) e - a{l+kc) ^ k)t . 

Let 



X l = r l+ '-T\ 



-i+l ^ 

SYi — I 

then 



k-l 



i=0 

where Xi are independent exponential random variables with rates a(l + ic)(n — i). By 
convolution, 

f T k+i(t) = / f T k(u)f T k+i_ T k(t — u)du 
t 

= a(l + kc)(n-k) ( f Tk (u)e- a{1+kc){n - k)(t - u) du, 
Jo 

for k — 1, 2, . . . , n — 1, where 

/ T i(i) = f T ^At) = nae~ nat . 

□ 
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Corollary 1 Assume that c ^ 1/i for % — 1, 2, . . . , n — 1. Then the unconditional density 
junction of r k is given by: 



k-1 



fr k (t) = J2 a K 



ae 



-fijCLt 



(7) 



3=0 



where the coefficients are given by: 

C*k,jPk 



Pk-Pj 

k-1 a 
LXk,uPk 

u=0 Pk ~ Pu 

Pj = (n-j)(l+jc) 



j = 0,l,...,A;-l 



j = k 



and ct\ o = n. 



Applying the recursive formula (4) iteratively gives the above corollary in which the 
analytic expression are given with recursive formulas for the coefficients. From the view 
point of computational convenience, we can see the advantage of the recursive formulas 
for the coefficients. The same result is also obtained in Zheng & Jiang (2009), with 



in _ kmk _ i _ j)!(n^ Wi (l + (m + j - n)c)) ' 

We remark that the method to derive the recursive formula of r fc stated here is related 
to the total hazard construction method adopted by Yu (2007) and Zheng & Jiang (2009). 
Assume we first enter the market right after the /cth default of the n entities. Then r k+1 — 
r k is the first default time being observed. We draw a collection of independent standard 
exponential random variables E±, E 2 , . . . , E n _^. By using the total hazard construction 
method, we have 



T k+i _ T k = min 



E 



n—k 



a(l + kc) ' a(l + kc) ' ' a(l + kc) 



which implies that 



n—k 



P(t 



fc+i 



T 



> *) = n p ( E i > a ( x + kc ) 1 ) = e 



-a(l+kc)(n-k)t 



i=l 



We note that the computational cost for the density of r k can grow up quickly, when k 
is getting large. The following propositions demonstrate the behavior of r k when k and 
n are large. In Propositions 2 and 3, we temporarily define r k (n) as kth default time in 
a portfolio of n names. 



Proposition 2 For a fixed k, when n — > oo, r k (n) — > almost surely. 
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Proof: For any given e > 0, by Markov's inequality, 

E[r k (n)} 2 



P{r\n) > e) < 

From Proposition 1, 



f 2 



^n)]=g^ ] = | \ (1 + 4_ z) <^ 
Var[r k (n)] = yVar[X t ] = g . J < 



i=o i=o K 1 + ^c)(n-z)) 2 a 2 (n-£;) 2 ' 

Since 

oo oo k2 cxd l 

then from Borel-Cantelli Lemma, we have with probability 1, for all large n, r k (n) < e. 
Hence T k {n) — > almost surely. □ 

Proposition 3 Let k — > oo (due to n> k, in this case, n oo), r k (n) — > almost surely 
and in particular, T n {n) — > almost surely as n — > oo. 

Proof: From Proposition 1, 

- fc / \i ^„ rv i ^ 1 ^ 1 2c + 2 .„ ln(l + cJfe). 

e\t (n)} = y E\Xi\ = y — — -- — r < y — — ^ — - < — - - — - i+- - 

[ v n ^ o(l + ic)(n - i) a(l + ic)(k - i) a(ck + l) K c 



i=0 



Indeed, 



y = y — - — \^— + ^—} 

fr^ a{\ + ic)(k - i) f^ a(ck + l) L k-i 1 + d J 

1 1 c 

= y — - — \— + —^- 1 

^ a(c£; + l) L i + 1 l + d J 
2c+2 " 1 



< 



< 



V- 

a{ck + 1) ^ 1 + ci 

2c + 2 /* 1 

a(cfc + 1) ^ io 1 + cx Z ' 
2c + 2 ln(l + cfc) 

a(cfc + l)^ c 



Similarly, 



fe-l fe-l 1 fO^_|_9"\2 oo | 



^ (a(l + ic){k - z)) 2 a 2 (cfc + l) 2 ^ (1 + *c) 2 



Since 

\2 



~ ~ (2c+2) 2 . ln(l + cfc). 2 

and 

- ~ (2c + 2) 2 ~ 1 

y Var[r k (n)} < y \ ' y — < oo, 



fc = l fe: 



-a2(cfc + l) 2 ^(l + ic) 5 



we have 

oo 

$>[r fc W 2 ]<°c. 



k=l 

-kl 



Using the same argument in the proof of Proposition 2, we can deduce that r (n) — > 
almost surely when k — > oo. j-j 

The above two propositions give us some insight about the infectious contagion, when 
the portfolio size becomes large, with entities inside the portfolio being associated by 
infection, the contagion becomes more intensive. 

3.1 Stochastic Intensity 

In this subsection, we extend the homogeneous contagion intensity process to the case 
that the constant intensity rate a is replaced by an "exogenous" stochastic process X, 
i.e., 

Proposition 4 Suppose there are n entities in our portfolio, where the contagion stochas- 
tic intensity process follows (8). Then the unconditional density functions of r k , k = 
1,2, ... ,n, given the realization of (X(s))q< s<oq , satisfy the following recursive formula: 

fr^\(x {s))o<s< Jt) = (n- fc)(l + kc)X(t) ft ^|(x (s )) 0<a<oo (^)e- ( "- fc)(1+fcc) ^ X(s)ds ^, 

(9) 

for k = 1, 2, . . . , n — 1, where 



fr^x {s))Q< _ 3< Jt)=nX(t)e~ n I>^ d \ 

Proof: We note that in this case, r fc+1 — r k depends on the kill default time r k , in the 
way that the ik + l)th default rate X h+l (t) follows: 

X k+1 (t) = X(t)(l + kc){n - k) for r fe <t< r k+1 . 

Then we have the following key relationship 

P(r k + 1 -T k >t\r k ,(X(s)) < s<oo ) = e -f^ k+ Hr k +s)ds 

_ e -(l+kc)(n-k)f ( ' ) X(T k +s)ds 

Therefore, the density function of r k+1 — r k given the realization of r k and (X(s))q< 



s<oo 



IS 



dP(r k+1 - r k >t\r k , (X(s)) 0<s<oo ] 

/t*+i-t*|t*,(x(«)) 0S . <oo ^; 



dt 



(n - k){\ + kc)X(r k + t) e - (n - k)il+kc) So x (- k + s ) ds 
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Thus the density function of r fc+1 given the realization of {X(s))q< s<00 , is given by 

fr^\(X(s)) < s<oo (t) 

= Jo fr*\(x(s)) 0<3< Ju)(n - k)(l + kc)X(t)e-^ 1+kc ^r x (^du (10) 
= (n-k)(l + kc)X(t) J* ^|(x(,)) , s<oo (^e" ( "- fc)(1+fcc) ^ XWd ^n, 

for k — 1, 2, . . . , n — 1, where 



fr W s ))o< _ s< Jt) = fr^rO l{X (s )h< _ a< Jt) = ^W^^^. 

□ 

Again, applying the recursive formula (10) iteratively gives the following corollary, 
which was obtained in Zheng & Jiang (2009). 

Corollary 2 Assume that c ^ 1/i for i = 1, 2, . . . ,n — 1. Then the density function of 
r k , given J 7 *, is given by 

fr*\(x (s))o<a< Jt) = J2a k>j X(t)e^I>^ du , (11) 

3=0 

where ctkj and f3j are given in Corollary 1. 



4 Multi-state Stochastic Intensity Process 

In this section, we consider the stochastic intensity process to be a multi-state Markov 
process. Here the state space S — {1,2, ... , N} represents the set of all exogenous states. 
For simplicity, we reduce the number of states to two. Then the stochastic intensity 
process X(t) alternates between x\ and Xi, so that 



X(t) 



X\, when the exogenous state lies in 1 
X2, when the exogenous state lies in 2. 

Let r)i denote the rate of leaving state i and 7Tj the random time to leave state i, where ni 
is an exponential random variable, i.e., 

p( n . > t) = e- Vit , i = 1,2. 

In this case, we consider a two-state, Markovian regime-switching, intensity-based model 
for portfolio default risk. Guo (1999) considers a two-state, Markovian regime-switching, 
model for option valuation. Here we follow her idea to derive the unconditional density 
function of r fc . Let Ti(t) be the total time between and t during which X(s) = x±, 
starting from X(0) = X\. We then draw exponential random variables £i with intensity 
rji and £2 with intensity ?7 2 independent of X. Therefore 

HWH 6 + T f ~ U ft* T 2 m{ T «-^ (12, 

t, £1 > t U, C,2 > t 
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where "=" means "equals in distribution". Let 

ipi(l,t) = E(e- im) ), i = 1,2, I eTZ. 

By Equation (12), we have 

V>i(7,t) = JoVie~ im+l)u M 1 ^ ~ u)du + e~^+ l )\ 
Ml, t) = Jo r/ 2 e-" 2 >i(/, t - u)du + e"" 2 *. 

Taking the Laplace transform on both sides gives: 

C(Mht))(; s) = + T^CiMhtm, s), 

£(^ 2 (Z,t))(., s ) = -J- + -^_£(^(M))(., S ). 

Therefore, 

rt i fi A\( \ s + r/! + r/ 2 



S 2 + S% + S7] 2 + Si + lr] 2 ' 

By taking the inverse Laplace transform of the above equation, we have 

e- Q '[cos(v^;t) + sm(y/ut)], to > 

e- Q '(l + /3t), u = 

e- a '[cosh(v/^cJt) + -^sinh( v / = a;t)], u < 0, 



where 



V1+V2 + I Q V1 + V2-I , 2 

a = , p = , u = lrj2-a . 



Similarly we have 

C{ifa{l,t)){;s) = — ■ , 

+ S7/i + S7/2 + S/ + £?72 

and 

{e- a '[cos(v^Jt) + ^ sin( v /cJt)], w > 

e- Q '(l + at), u = 

e- Qi [cosh(v^cJt) + -^sinh( v /= wt)], w < 0. 

We proceed with the derivation of the unconditional density function of r k . Given X(0) = 
E(e- L ti x{s)ds ) = E(e- L ^™ +X ^- T ^) = e- Lx *^{L{x x - x 2 ),t), L > . 

Then 

E(X(t)e Jo ) = -_ . 

Combining the results in Corollary 2, the unconditional density function of r k can be 
obtained: when X(0) = Xj, 

fAt) = ~h~w it • (13) 
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5 The Heterogeneous Case 



In this section, we present our method in obtaining the unconditional distribution of the 
fcth default time in the heterogeneous case of the interacting intensity default process. 
For simplicity of discussion, we consider a two-group case. The first group(Gi) consists 
n\ obligors and Aj(i) denotes the default rate of name % in G\ at time t, while the second 
group (G2) consists obligors and Aj(i) denotes the default rate of name % in G2 at time 
t. The interacting intensity process of the two-group case is assumed as follows: 



1 + h H l {r j <t} + cEl{f j <t} 



Xi(t) = a 

Xi(t) = a [ 1 + bJ2 l {r J <t} + cJ2 l {fj< 



(14) 



t} 



where tj and fj denote the default time of name j in G\ and G2, respectively, a, a are 
positive constants and b, c, b, c are nonnegative constants. Let N l be the number of defaults 
in G\ right after the ith default of our portfolio, where N° is assigned to be 0. 

Proposition 5 Suppose our portfolio has two groups of entities G\ and G2, where G\ 
and G2 consist of ni and ri2 obligors, respectively, and n = ri\ + n^- For each t > 0, let 
\i(t) and \i(t) denote the default rates of name i in G\ and G2 at time t, respectively. 
These default rates follow (14)- For each i = 1,2, ■■■,n, let N % denote the number of 
defaults in G\ right after the ith default of our portfolio. Then the recursive formula of 
the joint distribution of T k and N h is given by: 



/ r fc+i jvfc+i (t, m + 1) 



+(k, m [ f Tk M^ m y~ iCk ' m+L ' m){t ~ u)d ^ 



(15) 



where 



f T k N k(t,m) 



dP(r k > t, N k 



dt 



and 



(k,m — a{ni — m)[l + bm + c(k — m)}, 
(k,m = a(n 2 — (k — m))[l + bm + c(k - m)], 

f T i tN i{t, 1) = riiae-( nia+n2<i )', 
f T i,m(t,0) = n 2 ae- ( - nia+n2 ^ t , 
f T i tN i(t,m) = 0, m^0,l. 

Proof: Note that both r h+1 — r k and N h+1 are independent of r fc , but depend on N k . 



Conditioning on N h , r 



k+l 



T 



c and N h+1 are independent. The (k + l)th default rate 



X k+1 (t) that triggers r fc+1 , given N k , is given by 

A fc+1 (t) = C k ^ N k + ( kN k for 



r k < t < r k+1 . 
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Hence, we have the following relations: 

p(r k+1 — r k > t \ N k ) = e~^ k ' Nk+ ^ k - Nk ^ 



and 



On the other hand 



f rk+ i_ T k lNk (t) = (Ck,N k + Cfe.ArOe Kw+tw*. 



P(N k+1 - N k = 1 | N k ) = P(the (k + l)th default happens in G x \ N k 

Ck,N k 



and similarly, 



then the result follows 



Ck,N k + Ck,N k 

P( iV fc + 1 - N k = | N k ) = ^ . 

Ck,N k + Ck,N k 



xnen xne resuu ronows. 

Corollary 3 Assume that p\ m 7^ f or 1 < k < n, max{0, fc — n 2 } < m < min{/c,ni} 
and i = 0, . . . , A; — 1, j = 0, . . . , m. Then the joint density function of r k and N k is given 



by: 



k—l m 

f T k^k{t,m) = ^2 "fc,m,i,j 



"ft,,* 



(16) 



i=0 j=0 



where the coefficients are given by the recursive formula: 

"fc,m+l,i,jCfc,m+l &k,m,i,jCk,m • _ n -i i i • _ n ™ 

— - — — h — — — , I — U, 1, . . . , K — 1, J — U, ... , 777, 

Pfc,m+1 Pi J Pk,m Pi,j 

at,m+i, ij -a,m + i « = 0, 1, . . . , A; — 1, j = m + l 



"fc+l,m+l,i,j — < 



Pfc,m+1 ~ 
k—l m+1 /• 

u: fc,m+l,u,i; l 5fe,m+l 

u=0 i>=0 fik,m+l At, 
fe-1 m a 



u=0v=0 P\m Pu,' 



0. 



% = k, j = m + 1 

i — k,j — m 
otherwise 



Ci,j Ci,j 

and the boundry conditions are as follows: 

"1,0,0,0 = ra 2 a, "1,1,0,0 = nia, ai,i,o,i = 0, ai, m ,ij = 0,m ^ 0, 1, 



and 



"fe+l,0,i,0 — * 



Cfc,0"fc,0,i,0 
p\o — A,0 



fe-1 



Cfc,0"fc,0,u,0 



i = 0. ...,k-l 



i = k. 



u =o Pkfi — Pu,0 

As a result, the unconditional density function of r k is given by 

min{fc,ni} 

fr"(t) = fr k ,N k (t,m). 
m=max{0,fc— 712} 
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6 Evaluation of Basket CDS and Sensitivity Study 

Consider a kth to default basket CDS with maturity T. Assume Sk is the kth swap rate, 
and R is the recovery rate and r is the annualized riskless interest rate. The protection 
buyer A pays a periodic fee S^A; to the protection seller B at time t i: i = 1,2, . . . , N, 
where = to < t\ < . . . < tjy = T and Aj = ti — £j_i. If the fcth default happens in the 
interval [tj,tj + i\, A will also pay B the accrued default premium up to r k . On the other 
hand, if r k <T, B will pay A the loss occurred at r fc , that is, 1 — i?. Then the swap rate 
Sk is given by 



We remark that one can reproduce the results in Zheng & Jiang (2007) with our proposed 
methods in a reasonably high accuracy. 

In Table 1, we present the swap rates under the intensity-based default contagion 
model with exponential decay. We compute the swap rates with n = 2, k = 2, by selecting 
different values of the parameters. One can see that when we fix a, d, the swap rate 
decreases while c decreases. And when a, c are fixed, the swap rate falls down while d 
goes up. When c, d are fixed, the swap rate increases when a does. 

Table 2 presents the swap rates under the two-state, Markovian regime-switching, 
intensity-based default contagion model. Condition 1 is equivalent to the homogeneous 
case as X\ — X2- And we can see that if X\ < X2, swap rates increase when t]i increases, 
while swap rates decrease when t/ 2 increases. 

We present the swap rates under the heterogeneous intensity-based default contagion 
model in Table 3, where there are two groups of entities and the intensity process follows 
(14). We conduct our numerical experiments via selecting different value of the param- 
eters. As shown in Table 3, Conditions 1, 3 are equivalent to the homogeneous case for 
b = b = c = c. Condition 2 means that defaults in each group have a significant impact on 
the entities of the same group, but have a marginal impact on the other group. Condition 
4 shows the infection caused by entities in G\ is much more intense than those in G*2- 
We also make comparison on the analytic pricing approach (AP) presented previously 
and Monte Carlo approach (MC) in this case. The analytic approach takes less than one 
second to compute all the swap rates of one column with MATLAB on a computer with 
an Intel 3.2 GHz CPU, while the Monte Carlo approach takes more than 5 minutes to 
run 100,000 simulations. 

To study the sensitivities of the swap rates to a change in underlying parameters, we 
presents the derivatives representing the sensitivities in the homogeneous case with the 
intensity (5). Note that 



Sk 



(l-R)E(e-^ k l {T ^ T} ) 



(17) 



Eila E(A ie -rUl {Tk>u} + (r* - ^_i)e— k l 



{t»_l<T fc <ti}) 
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k=l 

k=2 

k=5 

k=1 



-k.1 
k=2 
k=5 

- k=lo 



Figure 1: The derivatives of swap rates S& with respect to a [left:AP, right:MC] (n—10, 
c=0.3) 



and 



where 



dfAt) 



k-l 



dc 



j=0 



a k,j1k,j + a k^k,j, 



0,1, 



k-l 



a 



k+l,j 



k-l 



H a k,ul'k, u + a 'k,ulk,u, j = k 



u=0 



IK 



Ik, 



Pk-Pj 



'" 3 dc ' 

and a l = 0, a^j and j3j are given in Corollary 1. Then combining (17), One can have 

dS k , n t s dS k 



9 J a) 



and k (c) 



da dc 
We present the derivatives 6k(a), 6k{c) by using the recursive formula(AP) above (Figure 
1, 2 (Left)), while we also use Monte Carlo method (MC) and difference quotient to find 
derivatives (Figure 1, 2 (Right)), where we select step size as 0.1 in MC for difference 
quotient and 100,000 simulations have been done. We remark that using the Monte 
Carlo method to compute the derivatives is very time consuming, and the results are 
unsatisfying, i.e., 9i(c) is constant zero by definition, while the results are quite fluctuating 
by Monte Carlo. 



7 Concluding Remarks 

In this paper we propose a simple recursive method to compute the fcth default time 
distribution under the interacting intensity default contagion model (1). We simplify the 
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Figure 2: The derivatives of swap rates Sk with respect to c [left:AP, right :MC] (n= 
a=0.l) 



Table 1: Basket CDS rates with exponential decay intensity-based default model (n 
2,k = 2,T = 3, A = 0.5, R = 0.5, r = 0.05) 





c 


0.2 


1 


5 


a 


d 




0.1 


0.001 


0.0134 


0.0211 


0.0479 




0.01 


0.0134 


0.0210 


0.0477 




0.1 


0.0132 


0.0203 


0.0459 




1 


0.0123 


0.0160 


0.0322 




10 


0.0115 


0.0120 


0.0147 




100 


0.0114 


0.0114 


0.0117 


1 


0.001 


0.3654 


0.4961 


0.7529 




0.01 


0.3651 


0.4955 


0.7526 




0.1 


0.3626 


0.4898 


0.7502 




1 


0.3464 


0.4390 


0.7184 




10 


0.3262 


0.3447 


0.4392 




100 


0.3222 


0.3242 


0.3342 



0.15 



0.1 





























'. r 
/■'. . 


\ 








k=2 

- k=5 

- k=10 






\ 




- \ ^ - 
1 >\ > 

V v 
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Table 2: Basket CDS rates in two-state stochastic intensity process case (n = 10, T = 
3, A = 0.5, R = 0.5, r = 0.05, c = 3, X(0) = X\, Condition 1: X\ = l,x 2 = 1, rji = 1, r/ 2 = 
1, Condition 2: xi = l,x 2 = 2, 771 = 1,772 = 1, Condition 3: X\ = l,x 2 = 2, 771 = 1,772 = 2, 
Condition 4: x\ = l,x 2 = 2,771 = 2,772 = 1) 



; 


Condition 1 


Condition 2 


Condition 3 


Condition 4 


1 


5.0242 


5.2507 


5.2409 


5.4575 


2 


3.9288 


4.1170 


4.1087 


4.2891 


3 


3.4456 


3.6184 


3.6106 


3.7766 


4 


3.1369 


3.3005 


3.2930 


3.4503 


5 


2.9035 


3.0605 


3.0532 


3.2043 


6 


2.7070 


2.8588 


2.8516 


2.9979 


7 


2.5270 


2.6743 


2.6672 


2.8093 


8 


2.3473 


2.4904 


2.4833 


2.6214 


9 


2.1459 


2.2847 


2.2775 


2.4114 


10 


1.8608 


1.9945 


1.9870 


2.1159 



Table 3: Basket CDS rates in heterogeneous case with 2 groups {ri\ = 5,7i 2 = 5, T = 
3, A = 0.5, i? = 0.5, r = 0.05, Condition 1: a = a = l,b = b = c = c = 3, Condition 2: 
a = a = l,b = c = 3,b = c = 0.3, Condition 3: a = a = l,b = b = c = c = 0.3, Condition 



4: a 


= a = 1, 6 = 


b = 3, c 


= c = 0.3) 












fc 


Condition 


1 


Condition 


2 


Condition 


3 


Condition 


4 




AP 


MC 


AP 


MC 


AP 


MC 


AP 


MC 


1 


5.0242 


5.0265 


5.0242 


5.0352 


5.0242 


5.0205 


5.0242 


5.0463 


2 


3.9288 


3.9352 


3.4752 


3.4692 


2.7073 


2.7167 


3.2065 


3.2167 


3 


3.4456 


3.4510 


2.8287 


2.8245 


1.9036 


1.9123 


2.5866 


2.5922 


4 


3.1369 


3.1417 


2.4246 


2.4209 


1.4799 


1.4860 


2.2543 


2.2567 


5 


2.9035 


2.9062 


2.1161 


2.1135 


1.2081 


1.2095 


2.0302 


2.0333 


6 


2.7070 


2.7068 


1.8376 


1.8366 


1.0112 


1.0116 


1.8554 


1.8549 


7 


2.5270 


2.5270 


1.6445 


1.6392 


0.8550 


0.8535 


1.7036 


1.7013 


8 


2.3473 


2.3477 


1.4821 


1.4757 


0.7203 


0.7205 


1.5582 


1.5545 


9 


2.1459 


2.1440 


1.3215 


1.3171 


0.5921 


0.5920 


1.4015 


1.3985 


10 


1.8608 


1.8625 


1.1169 


1.1096 


0.4451 


0.4448 


1.1889 


1.1851 



17 



problem in the homogeneous case with exponential decay (2) and with Markovian regime 
switching stochastic intensity (8). We further consider the problem in a two-group het- 
erogeneous case (14). We then present the numerical results for the basket CDS rates and 
sensitivity study using the proposed method. The main advantage of this method is that, 
by using the kth default rate, one can deduce the distributions of kth default times by 
recursive formulas with which one can easily compute CDS rates. Moreover, the proposed 
method can also be applied to the heterogeneous case to obtain the analytic expressions 
of CDS rates. Another key advantage is that, one can have the analytic formulas of the 
derivatives of the swap rates to the underlying parameters. The analytic formulas are 
fast and accurate while the Monte Carlo method is slow and inaccurate as the numerical 
experiment reveals. 

Acknowledgment: The authors would like to thank Prof. Mark H.A. Davis for his 
helpful discussions and suggestions. 
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